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We discuss the nature and origin of a rich complex of narrow absorption lines in the 
quasar J102325. 31+514251.0 at redshift 3.447. We measure nine C IV(A1548,1551) absorp- 
tion line systems with velocities from -1400 to -6200 km s~ x , and full widths at half minimum 
ranging from 16 to 350 km s -1 . We also detect other absorption lines in these systems, in- 
cluding H I, C III, N V, O VI, and Si IV. Lower ionisation lines are not present, indicating a 
generally high degree of ionisation in all nine systems. The total hydrogen column densities 
range from <10 17,2 to 10 19,1 cm~ 2 . The tight grouping of these lines in the quasar spectrum 
suggests that most or all of the absorbing regions are physically related. We examine several 
diagnostics to estimate more directly the location and origin of each absorber. Four of the 
systems can be attributed to a quasar-driven outflow based on line profiles that are smooth 
and broad compared to thermal line widths and to the typical absorption lines formed in in- 
tergalactic gas or galaxy halos. Several systems also have other indicators of a quasar outflow 
origin, including partial covering of the quasar emission source (e.g., in systems with speeds 
too high for a starburst-driven flow), O VI column densities above 10 15 cm~ 2 and an apparent 
line-lock in C IV (in two of the narrow profile systems). A search for line variability yielded 
null results, although with very poor constraints because the comparison spectra have much 
lower resolution. Altogether (but not including the tentative line-lock) there is direct evidence 
for 6 of the 9 systems forming in a quasar outflow. Consistent with a near-quasar origin, eight 
of the systems have metallicity values or lower limits in the range Z Js 1 — 8 Z P . The low- 
est velocity system, which has an ambiguous location based on the diagnostics mentioned 
above, also has the lowest metallicity, Z $C 0.3 Z , and might form in a non-outflow envi- 
ronment farther from the quasar. Overall, however, this complex of narrow absorption lines 
can be identified with a highly structured, multi-component outflow from the quasar. The high 
metallicities are similar to those derived for other quasars at similar redshifts and luminosities, 
and are consistent with evolution scenarios wherein quasars appear after the main episodes of 
star formation and metal enrichment in the host galaxies. 
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1 INTRODUCTION 

Quasars represent episodes of rapid supermassive black hole 
(SMBH) growth and probably a unique period in the early evo- 
lution of galaxies. They may directly follow a major merger 
jPerez-Gonzalez et alJl2008klHopkins et alJlJ OOS) or a big blowout 
of gas and dust. However, the nature of the relationship between 
SMBH growth and galaxy formation is not well understood. Feed- 
back from quasar outflows may play an important role in the evo- 
lution of this relationship. We are using narrow absorption lines 
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(NALs) in quasar spectra to study quasar outflows and environ- 
ments across a range of scales. 

NALs have full widths at half minimum (FWHMs) less than 
several hundred km s , and they appear in a variety of ultravio- 
let (UV) resonance transitions, including C IV A1548 and A1551, 
N V A 1239 and A 1243, C III A 977, Si IV A 1394 and A 1402 and 
Ly a A 1216, Ly / 3 A 1026, Lyr A973 Lyg A950, and Lye A 93 8 
jFoltzet all [198(1 1 Anderson et all 1 19871 ; [Hamann & SaQl2004h . 
The first step in using these absorption lines to study the quasar 
environment is to determine simply where the lines form. There 
are several possibilities, including unrelated (intervening) clouds 
and galaxies which happen to lie in the line-of-sight, nearby clus- 
ter galaxies and intrinsic clouds within the quasar host galaxy 
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and i t s extended environment i Ganguly et_alJ 12001 1_ Vestergaard 
20031; iTrump et al.1 120061: iNestor et al.1 l200sf ; Iwild et al.1 |2008J ; 
Gibson et all2008h . 



The statistical excess of C IV NALs at velocity shifts v > 
— 12,000 km s , where negative v indicates motion towards the 
observer, indicates that many of these absorb ers are directly re- 
lated to quasar enviro nments ( Wevmann et al.l 1 19791 : INestor et al.l 
l2008l;IWild et alj2008l) . The excess is largest at v > - 1000 km s" 1 , 
where roughly 80% of C IV systems with rest equivalent width 
REW(1548A) ^ 0.3 A have a quasar-related origin (Nestor et al.l 
2008). This intrinsic gas can form in quasar-driven outflows, 
starburst-driven outflows, merger remnants or ambient gas in the 
host halos, or in other galaxy halos in the same galaxy cluster as 
the quasar. At higher velocities, the excess can be attributed di- 
rectly to quasar -driven outflows. In the velocity range -1000 to - 
12,000 km s-'. lNestoreTall l teOOSh estimate that > 43% of C IV 
NALs with REW(1548A) ^ 0.3 A originate in a quasar outflow. 

NALs encompass a wealth of information about the basic 
properties of quasar outflows, which we use to gain insights into the 
outflow physics, acceleration mechanisms and geometry of the near 
quasar environment. NALs represent a very different type of quasar 
outflow compared to the well-studied much broader and higher ve- 
locity broad absorption lines (BALs), but they might simply be dif- 
ferent manifest ations of a single outflow phenomeno n viewed at 
different angles JGangulv et all 1 999Ll200ll:lElvisl200Ch . The NALs 
that do not form in outflows probe the gaseous environments of 
quasars more generally, and can be used to examine conditions in 
different host galaxies environments such as starburst-driven out- 
flows or larger-scale gas distributions produced by messy mergers. 

We are involved in a program to study the location, origin and 
abundance information for absorbers in a sample of high redshift 
quasars. We are particularly interested in high redshift quasars be- 
cause z ~ 2 — 4 is the cosmic era when host galaxies are thought to 
grow ra pidly and form most of their stars, possibly throug h merger 
events jPerez-Gonzalez et al.l [20081 ; iHopkins et al.ll2008l) . Choos- 
ing redshifts above z ~ 2.7-3 also allows us to measure lines at 
shorter rest-frame wavelengths with ground-based telescopes, in- 
cluding importantly, the H I Lyman series, in order to obtain more 
and better constraints on the absorber ionisations, column densi- 
ties and metal abundance. We are interested in using abundances 
to discern rough star formation histories of quasar host galaxies in 
order to make inferences about the relationship between the quasar, 
the growth of the central black hole, and the evolution of the host 
galaxy. 

Broad emission lines (BELs) have been used most of- 
ten to study quasar abundances. The most reliable results sug- 
gest metallicities of at least solar, and up to a few times so- 
lar, which requires significant previous star formation in the 



violet (UV) spectral range, where many useful metal and Hydrogen 
lines occur. 

NALs offer certain advantages in the study of metallicities and 
other gas characteristics in the near-quasar environment. Their nar- 
row widths mean the C IV doublets, separated by 500 km s , are 
resolved. We use resolved absorption line doublets to disentangle 
saturation effects, and to obtain accurate line optical depth and col- 
umn density measurements. NALs also form in a range of physical 
locations, providing a more complete picture of the regions near 
quasars. Because the NAL methods are completely independent 
from the BEL methods, requiring only column densities and ion- 
isation corrections, the NAL metallicities provide an independent 
test of the BEL results. 

Here we present results for the particular luminous quasar 
J102325. 31+514251.0 (hereafter J1023+5142) at a redshift of 
Zem = 3.447, which is during the peak in quasar activity and the 
epoch of rapid galaxy formation. This quasar contains a rich com- 
plex of nine distinct C IV NAL systems at velocities from -1400 to 
-6200 km s _1 . The density and diversity of lines in this complex 
merits special attention. We will argue below that most (or all) of 
these systems form in a highly structured quasar-driven outflow. 

To interprete the metallicities and other data provided by these 
NALs, we examine several diagnostics that can identify intrinsic 
NALs that form i n quasar-driven outflows dHamann et alj|l997l : 
iBarlow et al.ll 19971) . In particular, 1) variability studies have found 
intrinsic absorbers varying on relatively short timescales of months 
to years, providing strong evidence for these absorbers belonging 
to outflows either crossing the line of sight to the quasar or expe- 
riencing changing ionisa tion with the variations in the continuum 
emission dHamann et al. 19971; Barlow et al. 19971 : lAldcroft et al.l 
ll997l:lNaravanan et alj|2004 iMisawa et alj|2007|) . 2) Detection of 
partial coverage of the background light source along the line of 
sight strongly implies gas forming very near the source. This phe- 
nomenon occurs when the absorbing 'clouds' are smaller than the 
background source, allowing part of the light from the source to 
reach the observer unabsorbed. This partial covering is easily de- 
tected in multiplets like the C IV doublet where the optical depth ra- 
tio between the two lines is fixed by the oscillator strengths. When 
the source is partially covered, some light fills in the bottom of 
the absorption line, and makes the apparent optical depth ratio ap- 
pear different than the real optical depth ratio. 3) Outflow lines 
tend to have profiles that are broad and smooth compared to ther- 
mal widths dHamann & Ferland! 1 19991; I Srianand & Petitiearfood: 



lGangulvetal.ll2006l : ISchaveetal.ll2007h . In well studied NALs. 
these three indicators (variability, partial covering and broad pro- 
files) tend to appear together, which further increases the proba- 
bility that the occurrence of an individual indicator accurately pre- 
dicts an outfl ow origin very near the q uasar for a given absorption 
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2003J ; IWarner et alj |2004| ; iNagao et alJ 120061 ; ISimon & Hamanri 
20ld) . The metal-ric h BEL result is true even for the highes t 
redshifts studied, e .g. IPentericci eTail d2002h : Ijiang et al.1 J2007h : 
ljuarez et al.1 Q2009J), with redshifts out to 6.4. The most reliable 
results based on BAL column de nsities suggest m etallicity ranges 
between solar and ten times solar dArav et alj200ll) . Previous stud- 
ies of NALs in low redshift samples have found super-solar metal- 
licities and highly ionised gas, and have successfully probed sev- 
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2004; Ganguly et al. 2006; Gabel et al. 2006). These lower redshift 



samples cover a wide range of luminosities, observed in the ultra- 



also note that super-solar metallicities are consistent with an in- 
trinsic origin for the gas. There are examples of high-metallicity 

as in intervening systems, but not of low- metallicity intrinsic gas 

Prochaska et alj2006l : ISchave et alj|2007l) . 

We describe the data acquisition and reduction in §|2] the iden- 
tification and fitting of the absorption lines in § |3.1| and |3.2| and the 
abundance and ionisation analysis in § 13.31 We briefly describe in- 
dividual absorption line systems in § [4] We discuss the arguments 
for the locations, probable intrinsic origins and quasar-driven out- 
flow properties of the gas in §[5] and conclude with a summary in 

m 
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Figure 1. Region of the spectrum of J1023+5142 with C IV absorption. Individual C IV doublets are labelled by number. The lower x-axis is observed 
wavelength in Angstroms, while the upper x-axis is velocity shift of the shorter wavelength doublet line at 1548.20 A from the quasar rest frame in kilometres 
per second. The flux units are normalised so that the continuum has a value of one. The gap between 6788 and 6797 A is a gap between Echelle orders in the 
spectrograph. The longer wavelength line of system 8 falls in this gap. Strongly blended lines are considered components of a single system, e.g. systems 7 
and 9. 



2 OBSERVATIONS AND DATA REDUCTION 

We observed the quasar J1023+5142 on March 29, 2007 with 
the Keck I HIRESr Echelle spectrograph as part of an observing 
campaign to measure spectra of high redshift quasars with known 
narrow associated absorption lines. We used an 0".86 wide slit 
for a spectral resolution of R ~ 40,000 or velocity resolution of 
~ 7 km s . Our data span the wavelength range from 3700 to 
8100 A corresponding to 830 to 1820 A in the quasar rest frame. 
This spectral range covers a variety of interesting lines, including 
rest-frame Ly y 970 A, C IV 1548, 1551 A and several other lines in 
the H I Lyman series down to the Lyman limit at 912 A. We use four 
exposures totalling 2 hours on the source. The spectral region from 
~ 3700 to ~ 4980 A is well-covered with considerable overlap be- 
tween Echelle orders at some wavelengths, but above 6540 A there 
are small 15-40 A gaps between orders (one such gap is apparent 
in FigureQ}, and two larger gaps at 4980-5070 A and 6575-6670 A 
where the spectrum falls into a physical gap between detectors. 

We reduce all the data using the MAKEE HIRES data re- 
duction package. The spectra are sky background-subtracted and 
extracted from the 2D frame with a low order polynomial trace us- 
ing a white dwarf standard observed the previous night with similar 
seeing conditions. We use a Thorium Argon (ThAr) lamp spectrum 
for wavelength calibration. The resulting spectra are on a vacuum 
and heliocentric wavelength scale. The spectra are not absolute flux 
calibrated. 

We normalise the spectra to unity by fitting a pseudo- 
continuum to all of the quasar emission, including the emission 
lines. The pseudo-continuum is defined as follows: for regions with 
few absorption lines, we apply a polynomial fit to the local con- 
tinuum in each Echelle order. We accomplish the continuum nor- 
malisation in crowded regions where the continuum is affected by 
significant absorption, e.g. the Ly a forest, by first averaging to- 
gether several adjacent spectral orders into a single spectrum. Then, 
we visually inspect the region for small sections of continuum not 
affected by absorption or obvious noise spikes, and interpolate be- 
tween these sections, fitting the entire region with a low order poly- 




4400 



4500 



4600 
Wavelength (A) 



4700 



4800 



Figure 2. Region of Ly a forest spectrum with the continuum fit over- 
plotted. The regions spanning the Ly j3 and O VI NALs are labelled above 
the spectrum. 



nomial. Our continuum fit for a region of the Ly a forest containing 
the Ly/3 and O VI NALs is shown in Figure[2] The continuum place- 
ment has an uncertainty of ~ 10% in the forest and 2-3% at other 
wavelengths. 



3 ANALYSIS 

3.1 Identification 

The broad, flat shape of the emission features in the spectrum of 
J1023+5142 make an accurate emission redshift difficult to deter- 
mine. The redshift provided by the SDSS spectrum is z em = 3.447. 
We estimate the reliability of this value by measuring the redshifts 
of the C IV A 1549, C III] A 1909 and Si IV+O IV] 11398 emis- 
sion lines using measurements of their centroids. We shift each 
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centroid respectively by -824, -730 and +36 km s to correct for 
known offsets from the nominal quasar redshi ft ([O III] A 5007 
emission), based on measurements by Tytler & Fanl d 19921) and 
IShen et alj (2007J) for average quasars, where negative values are 
blueshifts. The average redshift obtained from these emission lines 
is Zem = 3.429, which is offset from the SDSS value by Az = 0.018 
or ~ —1200 km s . These results and our efforts to measure 
the inherently uncertain emission line centroids suggest that the 
uncertainty of the redshift measured by SDSS is not more than 
Az em ^ 0.02, corresponding to Av ^ —1350 km s . We adopt the 
SDSS value throughout the remainder of this paper. 

We identify nine distinct C IV absorption line systems within 
6200 km s of the quasar redshift. We will refer to these as sys- 
tems 1-9, as indicated in Figure [TJ and Table [T] below. Other C IV 
systems are present at -16,800 and -33,800 km s _1 in the spec- 
trum, but they have narrow widths, complete covering, and blend- 
ing problems in the Ly a forest, which, along with their high veloc- 
ities, make them likely candidates for intervening gas and exclude 
them from further analysis in this work. 

After identifying the C IV doublets, we search the spectrum 
for other common NALs such as Si IV, N V, C III, O VI, and H I 
Lyman series lines at the same redshift. We also search for lower 
ionisation species, such as C II and Si II, but find none. All of 
the systems, except possibly system 1, appear to have relatively 
high ionisations based on the presence and absence of high and low 
ionisation species respectively. Each set of absorption lines at one 
redshift is considered a system, as labelled in FigureQ] Several of 
these systems are blends of two or three components, which are not 
individually labelled in the figure. 

Systems 1 and 2 (Figures[3]and[4]l appear to be line-locked in 
C IV. The velocity offset between the A 1548 line in system 1 and 
the A 1551 line in system 2 is remarkably small (< 2 km s _1 ) com- 
pared to the FWHMs of these lines (~ 30 km s _1 ) and the velocity 
shifts from the quasar systemic, ~ —1440 and ~ —1940 km s . 
If this overlap between the C IV lines in systems 1 and 2 repre- 
sents a physical line-lock, where the velocities of the two systems 
are actually separated by exactly their d oublet separation, and not 
a chance alignment in t he spectrum (see lGangulv et alj ( 120031) and 
iBraun & Milgroml dl989l) for full discussions of the possible physi- 
cal nature of line-locking), then it provides evidence for these lines 
forming in a quasar outflow driven by radiation pressure (see § 15.11 
below). 

3.2 Line Fitting 

We fit each NAL system with a Gaussian optical depth profile. The 
narrowest absorption lines are at least 1.5 times broader than the 
spectral resolution and the other lines are significantly broader than 
this. The absorption lines are, therefore, fully resolved, and such 
Gaussian optical depth profile fits are sufficient to determine accu- 
rate optical depths and covering fractions. The optical depths and 
covering fractions are held constant across the width of each line 
profile. Gaussian fits are actually essential to distinguish individ- 
ual absorption features in the crowded Ly a forest, and also useful 
to disentangle blended absorption in other areas of the spectrum. 
Furthermore, Gaussian fits smooth over noise spikes and large op- 
tical depth and covering fraction uncertainties in the wings of the 
lines. We also use Gaussian fits to simultaneously fit and lock to- 
gether various parameters including redshift, doppler b parameter, 
covering fraction and a 2: 1 optical depth ratio based on oscillator 
strength ratios for doublets such as the C IV, Si IV, N V, and O VI. 
To measure accurate optical depths, we must consider the pos- 



sible effects of partial coverage of the emission source by the ab- 
sorbing gas. The line of sight covering fraction affects the observed 
line intensity as follows: 



h = (l-C f )I +C f I e 



-T V 



(1) 



where ^ Cf ^ 1 is the velocity dependant line of sight cover- 
ing fraction, Iq is the emitted (unabsorbed) intensity and /y and Ty 
are the observed intensity and line optical depth at each velocity 
shift v. This equation assumes that the background light source has 
a uniform brightness given by In and the foreground absorber is 
homogeneous with a sing le value of Ty. T he viability of this as- 
sumpt ion is discussed by Haman n & Sabrai (2004) and lArav et alj 
(2005). We assume that all lines in a given multiplet have the same 
C/ at a given velocity. We do not explicitly attempt to distinguish 
between partial cover ing of the continuum source and of the BEL 
region as discussed bv lGangulv et alj Jl999D . However, we estimate 
from the SDSS spectrum that the C IV BEL peaks 20% above the 
continuum, which implies that the BEL can only account for partial 
covering of 0.8 or higher. 

We attempt to fit each system with the smallest possible num- 
ber of Gaussian components. This minimizes the number of free 
parameters and provides a more robust characterization of column 
densities, ionisations and abundances in absorbing regions whose 
internal velocities might be more complex than simple Gaussians 
(e.g., in outflows). We fit each absorption line with a single Gaus- 
sian unless 1) the system clearly has multiple components distin- 
guished by inflection points that stand out significantly above the 
noise fluctuations in the spectrum (e.g. system 7), or 2) a single 
Gaussian would miss a significant fraction, > 25%, of the absorp- 
tion line strength (e.g. system 4). The ~ 25% threshold is somewhat 
arbitrary, but it ensures that we achieve a good fit to the observed 
line and that significant portions of absorption (i.e., large enough 
to change the column density measurements) are not missed. For 
these exceptional cases we use the minimum number of Gaussians 
possible to achieve an accurate fit to the data. If a system is fit with 
two or more Gaussians, each Gaussian is labelled as a component. 
We assume that the covering fraction is the same for all components 
in a given system, such that the optical depths in Equation[T]simply 
add together in regions of component overlap (see Hamann et al. in 
preparation for further discussion). This simplifying assumption is 
well justified by the excellent fits to all the systems, with the possi- 
ble exception of system 6, which we discuss in more detail in §|4] 
All ions with Gaussian fits are shown along with their Gaussian 
optical depth profiles in Figures l3l9l Badly blended members of, 
e.g., the Lyman-series lines are not used to constrain ionisation or 
abundance in §[33] and are not shown. We check the Gaussian re- 
sults using direct integration and point-by-point measurements of 
Ty and Cf across the line profiles for several systems with either 
non-Gaussian profiles (System 8) or Cf < 1 (Systems 5 and 6) as 
discussed later in this section. 

The free parameters in the Gaussian fits are Cf, central optical 
depth, t, the Doppler b parameter, and redshift for each compo- 
nent in each system. We fit the C IV absorption lines first, then 
base the fit for other absorption lines on the C IV fit parameters. 
To ensure we are analysing the same gas in different ions, we fix 
the redshift for all absorption lines in a system to the C IV redshift. 
We further exclude unwanted contributions from blends or complex 
multi-phase gas by fixing the b parameter of all ions with ionisa- 
tions lower than C IV to that of C IV. Higher ionisation lines, such 
as N V and O VI, are allowed to be broader. However, we cap the 
b- values of the H I profiles at 140% of C IV. This cap is important 
for the abundance analysis below (§[33} because it ensures that the 



Narrow Absorption Line Complex at z = 3.45 5 



derived H I column densities do not include gas with dramatically 
different kinematics than C IV. We choose to cap the H I Z?-values 
at 140% of the C IV i-values instead of the much higher percent- 
age expected for purely thermal broadening because the widths of 
the C IV lines exceed the thermal widths expected for a gas pho- 
toionised by either the quasar or the inter-galactic UV spectrum. 
Therefore, we assume the fc-values are dominated by non-thermal 
broadening effects. On the other hand, setting the cap at 140% in- 
stead of something smaller, such as 100%, allows for some contri- 
bution of thermal broadening to b in the narrower systems (which 
would affect H I more than C IV). Overall, our fits to the Lyman 
lines should lead to reasonable but generously large estimates of 
the amount of H I gas that coexists with C IV, and therefore, to 
conservatively low estimates of the C/H abundance. 

As stated above, the covering fraction is a free parameter in 
the Gaussian optical depth fits of each doublet. In cases where the 
best fit profile has Cf < 1, we repeat the fit with Cf = 1 to test 
the robustness of the Cf < 1 result. We then compare how well 
each of the two fits with different values of Cf match the data. In 
cases where the fits are comparable, we assume Cf = 1, otherwise, 
the best fit is chosen. For example, we confirm that the Cf < 1 fit 
follows the data in systems 5 and 6, as shown in Figure(6] whereas 
the Cf = 1 fit does not match the observed doublet ratios in C IV 
and N V. We assume Cf = 1 for all singlet lines. The covering 
fraction for the Lyman lines is fixed at the C IV doublet covering 
fraction. This is necessary because the observed line ratios within 
the Lyman series are too severely affected by blending in the Ly a 
forest to yield their own independent measures of Cf. 

Table [7] lists fit parameters for all of the absorption lines that 
yield useful constraints for the ionisation and abundance analysis 
described in § 13.31 Absorption lines that are badly blended are not 
used in subsequent analysis, and are not listed in the table. Each 
system is listed separately by redshift and velocity shift, where neg- 
ative velocities denote gas moving towards the observer. Only the 
stronger member of the O VI doublet is listed, as O VI is never 
used in the abundance analysis because of either line saturation or 
strong blending in these lines in all systems. However, the strength 
of O VI is still useful as an indicator of the ionisation of the gas. 
In systems where N V is not present, we list upper limits for the 
stronger member of the N V doublet for completeness. TableQ]lists 
the central wavelength (X) and doppler b parameter values along 
with column densities and rest equivalent widths (REW) derived 
from the Gaussian optical depth profile fits. Systems 4 and 7 each 
have two blended components. The values of X,b, and logW are 
listed separately for these components in Table Q] but the REWs, 
listed only with the first component data, apply to the entire blend. 
We measure upper limits on H I column densities in all cases where 
all the Lyman series absorption lines are blended with intervening 
absorption lines in the Ly a forest. The same is true for singlet ions 
with upper limits on the column densities. 

We estimate uncertainties for the column densities by placing 
the continuum at the top and bottom of the noise around the fitted 
continuum, corresponding to the reasonable maximum/minimum 
values (~ 3 <J uncertainties) for continuum placement. We mea- 
sure 3 cr uncertainties for the H I column densities of 0.18 dex 
on average. The covering fraction is Cf = 1 for all the systems, 
unless otherwise noted in the footnotes of Table [T] e.g. systems 5 
and 6. We estimate the uncertainty in the covering fraction derived 
from the Gaussian profile fits both formally and informally. The 
formal uncertainties are estimated by propagating the error spec- 
trum through the calculation of Cf. However, these uncertainties 
are much smaller than the informal uncertainties, which are domi- 
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-2120 


Ly 7973 


4297.29 


0.035 


13.9 


s: 14.26 




Ly a 1216 


5371.61 


0.186 


13.9 


s: 14.26 




N V 1239 


5473.90 


;£0.006 


9.5 


s; 12.45 




Si IV 1394 


6158.52 


0.023 


7.2 


12.48 




Si IV 1403 


6198.30 


0.013 


7.2 


12.48 




C IV 1548 


6840.94 


0.096 


9.5 


13.62 




CIV 1551 


6852.29 


0.062 


9.5 


13.62 


4 3.41775 


Ly e 938 


4142.98 


0.059 


15.4 


sj 14.86 


-2182 


Ly 7973 


4296.43 


0.164 


15.4 


sj 14.86 




Ly a 1216 


5370.53 


0.589 


15.4 


s? 14.86 




NV1239 


5472.79 


<0.055 


11.0 


sj 12.28 




Si IV 1394 6157.28 


0.081 


9.1 


12.76 




Si IV 1403 


6197.09 


0.046 


9.1 


12.76 




C IV 1548 


6839.56 


0.215 


11.0 


13.56 




CIV 1551 


6850.91 


0.132 


11.0 


13.56 


3.41747 


Ly £ 938 


4142.72 


* 


46.8 


s; 14.64 


-2200 


Ly 7973 


4296.15 


* 


46.8 


sC 14.64 




Ly a 1216 


5370.19 


* 


46.8 


SC14.64 




NV1239 


5472.44 


* 


33.4 


^13.45 




Si IV 1394 6156.89 


* 


14.0 


12.75 




Si IV 1403 


6196.66 


* 


14.0 


12.75 




C IV 1548 


6839.12 


* 


33.4 


13.67 




CIV 1551 


6850.47 


* 


33.4 


13.67 



nated by uncertainty in continuum placement. We estimate cover- 
ing fraction uncertainties due to continuum placement uncertain- 
ties by first shifting the continuum near each red doublet mem- 
ber up and down by the 3rr continuum uncertainty, and then fitting 
the doublets with this new continuum. The actual uncertainties are 
probably smaller than the uncertainties we derive in this way, be- 
cause a similar shift in the continuum around both doublet members 
(a more likely occurrence) produces smaller changes in Cf. We find 
Cf = 0.7 ±0.15 for system 5 and Cf = 0.7 ±0.20 for system 6. If 
we fix the covering fraction in systems 5 and 6 at Cf = 1 instead of 
at the measured values, the column densities in all ions decrease by 
an average of 0.25 dex. 

We perform a simple test to determine the reliability of the 
Cf < 1 result from the Gaussian profile fits for systems 5 and 6. We 
predict the shape of the longer wavelength C IV and N V doublet 
members based on the intensity of the shorter wavelength member, 
combined with the 2: 1 T-ratio derived from the oscillator strengths 
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TableQ] Continued 












5" 3.40391 


C III 977 


4302.71 


0.059 


19.7 


=£13.05 


-3121 


O VI 1032 4544.53 


0.275 


25.3 


=£14.99 




Ly a 1216 


5353.70 


0.176 


32.2 


13.91 




NV 1239 


5455.65 


0.098 


13.0 


14.30 




NV 1243 


5473.18 


0.075 


13.0 


14.30 




CIV 1548 


6818.13 


0.141 


15.0 


14.12 




CIV 1551 


6829.45 


0.105 


15.0 


14.12 


6* 3.40196 


Ly y 973 


4281.07 


0.055 


31.7 


14.41 


-3254 


Ly ft 1026 


4515.18 


0.137 


31.7 


14.41 




OVI 1032 4542.51 


0.739 


63.5 


=£15.54 




Ly a 1216 


5351.33 


0.382 


31.7 


14.41 




NV 1239 


5453.24 


0.284 


50.0 


14.57 




NV 1243 


5470.76 


0.188 


50.0 


14.57 




CIV 1548 


6815.11 


0.220 


45.0 


14.02 




CIV 1551 


6826.42 


0.128 


45.0 


14.02 


7 3.39936 


Ly S 950 


4178.26 


0.008 


27.3 


=£13.19 


-3430 


C III 977 


4298.26 


0.102 


19.8 


=£12.36 




N III 990 


4354.49 


0.023 


19.8 


=£ 13.04 




Ly a 1216 


5348.17 


0.295 


27.3 


=£13.19 




NV 1239 


5450.02 


0.310 


18.0 


13.50 




NV 1243 


5467.53 


0.180 


18.0 


13.50 




CIV 1548 


6811.09 


0.187 


19.9 


13.09 




CIV 1551 


6822.40 


0.101 


19.9 


13.09 


3.39835 


Ly 8 950 


4177.31 


* 


60.0 


=£13.74 


-3496 


C III 977 


4297.28 


* 


43.7 


=£13.19 




N III 990 


4353.49 


* 


43.7 


=£13.04 




Ly P 1026 


4511.48 


* 


60.0 


=£13.74 




Ly a 1216 


5346.95 


* 


60.0 


=£13.74 




NV 1239 


5448.77 


* 


52.4 


14.22 




NV 1243 


5466.27 


* 


52.4 


14.22 




C IV 1548 


6809.53 


* 


43.7 


13.62 




CIV 1551 


6820.83 


* 


43.7 


13.61 


8 3.37983 


Ly p 1026 


4492.48 


0.135 


149.0 


=£14.29 


-4763 


VI 1032 


4519.68 


1.705 


149.8 


15.86 




Ly a 1216 


5324.43 


0.766 


149.0 


=£14.29 




NV 1239 


5425.82 


1.045 


155.5 


14.88 




NV1243 


5443.26 


0.643 


155.5 


14.88 




C IV 1548 


6780.85 


0.807 


155.5 


14.41 


9 3.35922 


Ly y 973 


4239.50 


0.039 


58.4 


=£14.02 


-6186 


C III 977 


4258.94 


=£0.158 


50.0 


=£13.35 




Ly P 1026 


4471.34 


0.112 


58.4 


=£14.02 




O VI 1032 4498.41 


0.413 


42.2 


14.51 




Ly a 1216 


5299.37 


0.513 


58.4 


=£14.02 




NV 1239 


5400.29 


0.212 


40.0 


13.84 




N V 1243 


5417.64 


0.115 


40.0 


13.84 




C IV 1548 


6748.94 


0.156 


40.4 


13.35 




CIV 1551 


6760.14 


0.083 


40.4 


13.35 



" System 5 has covering fractions C f (H i) = 0.70, C f (N v) = 0.67 
and C/(C IV) = 0.70. 

b System 6 has covering fractions C/(H i) = 1.0, C/(iVv) = 0.67 
and C/(C IV) = 0.72. 



elude that Cf < 1 for at least some portion of each of these lines in 
systems 5 and 6. 

The Gaussian fitting technique described above assumes a sin- 
gle covering fraction across the entire line. However, we demon- 
strate with the above T-ratio analysis that covering fractions are not 
always constant across a single line profile. To better account for 
this, and to determine if the guassian fits find reasonable average 
values for Cf, particularly for lines with very non-Gaussian shapes, 
we use a point-by-point method in addition to the Gaussian fitting 
method to determine Ty and Cf across the line profiles in several 
systems. We fit systems 5 and 6, the two systems with Cf < 1, 
along with system 8, which has a very non-Gaussian profile. We 
step across the absorption line, calculating average intensity in each 
small (a few times the resolution) regularly spaced sections of the 
spectrum, using the ratio of the intensities in the doublet in Equa- 
tion Q] to measure Cf and T v at each step. The point-by-point fits 
are shown in Figures QT] and [T2] The solid curve shows the shorter 
wavelength doublet member, while the dot-dashed curve shows the 
longer wavelength doublet member. The covering fraction at each 
point is represented by the filled circles. The steps used for sys- 
tem 6 are three resolution elements wide, which is wide enough 
to smooth over the noise but narrow enough to avoid blending the 
wings and core of the line. However, system 5 is narrow enough 
that using bins three resolution elements wide, or wider, across the 
wings of the line would blend too much information from the core 
and the continuum. Furthermore, the spectrograph resolution could 
be blending the covering fraction in the wings of the line with the 
continuum. Thus, for the narrow system 5, we measure only the 3 
resolution element bin at the line core. The step size for system 8 
is four resolution elements. This larger step size further smooths 
over noise, and can be used because the line is much broader than 
the other systems, lessening the impact of blending of the core and 
wings of the line. We derive formal covering fraction uncertain- 
ties (o"c,X represented as error bars at each point in FiguresfTTIand 
[72l The average C IV and N V central covering fraction from the 
point-by-point method matches the C IV and N V covering fraction 
derived from the Gaussian fitting method to within 10% in systems 
5 and 6. 

Based on this result, the results of the f-derived doublet ra- 
tio analysis, and uncertainties derived from the Gaussian fitting 
method, we are confident of the accuracy of the Cf < 1 measure- 
ments for systems 5 and 6, although the exact value of the covering 
fraction remains unknown beyond the fact that it is below 1. The 
C IV and N V column densities found in system 6 using the point- 
by-point method match the C IV and N V column densities derived 
from the Gaussian fitting method to within 0.14 dex and 0.06 dex 
respectively. The same comparison for N V in system 8 yields a 
difference of 0. 1 8 dex between the two methods. We conclude that 
the Gaussian technique is sufficient for comparing different sys- 
tems and generally provides accurate column density and covering 
fraction results for the purposes of this work. 



of each line. The predicted shorter wavelength member will only 
match the data if Cf = 1 . These predictions are shown in FigurefTol 
The observed data for both doublet members are plotted with solid 
curves and the predicted longer wavelength doublet member is plot- 
ted with a dot-dashed curve. The predicted shape of the longer 
wavelength member of C IV and especially of N V is much weaker 
than the observed shape for system 5. In system 6, the line centres 
of C IV, and more clearly N V, are stronger in the observed data for 
the longer wavelength members than in the predictions. We con- 



3.3 Ionisation and Abundances 

The abundance ratios can be derived from the ratio of measured col- 
umn densities corrected for the degree of ionisation in the gas. For 
example, the relative carbon to hydrogen abundance normalised to 
solar is given by: 



log 



N{C IV) 
N(HI) 



, ( /(go 



(2) 







Narrow Absorption Line Complex at z = 3.45 



where / is the ionisation fraction of a given ion, N is the column 
density and the final term on the right-hand side is the logarithmic 
solar a bundance ratio of hydrogen to carbon listed in lGrevesse et al.l 
(2007). The second term on the right is the ionisation correc- 
tion (IC). These correction factors can be large when comparing 
a highly ionised metal like C IV to H I. The exact values depend 
on the ionisation mechanism. Photoionisation by the quasar spec- 
trum is by far the most likely scenario based on the arguments in 
§ |5.1| that all of the systems are likely to be intrinsic to the quasar 
environment. 

We d erive values of IC using the photoionisation calculations 
shown in Haman n et al.l (2010). Their calculations adopt a nomi- 
nal quasar spectrum consistent with recent observational estimates 
at the critical ionising (far-UV) photon energieo The calculations 
also assume that the absorbing gas is optically thin in the Lyman 
continuum, which is appropriate for the column densities we mea- 
sure in the absorption lines of J1023+5142 (Table[T|(. 

Ideally, we would constrain the absorber ionisations by com- 
paring the ratios of observed column densities in different ions of 
the same element, such as N(C \u)/N(C iv) or N (N Ul)/N(N v), 
to the theoretical results in Hamann et al.l yOlOj). However, these 
constraints are only marginally usable in our data because N(C III) 
and N(N III) are always blended in the Ly a forest and are there- 
fore only ever constrained as upper limits. Therefore we estimate 
the IC from ratios such as N(N v)/N(C iv) or N(Si IV)/JV(C iv), 
with the additional assumption that the relative metal abundances 
are approximately solar. The specific ionisation constraints used for 
each system sometimes lead to upper limits, lower limits or specific 
values for the abundance ratios, and are described in more detail 
for individual systems in §|4]below. Our best estimates for the C/H 
abundances based on these constraints are all super-solar, except in 
system 1. Table|2] which contains several different abundance indi- 
cators for each NAL system, lists these estimated 'best' abundances 
in column 3, titled [C/H]i, est , for the nine systems. 

We also calculate robust lower limits on the metal to hydro- 
gen abundance ratios by applying minimum values of the ionisa- 
tion correction (ICmm . lHamann et al.l d 1997b ) to the measured C IV, 
Si IV and N V column densities, when available. Each metal ion 
has a unique global IC m ,„ that occurs near the peak of its own 
ionisation fraction. For example, f(H l)/f(C iv) peaks approxi- 
m ately where f(C iv) is largest. We use the values of IC„„ n listed 
in lHamann et alj J2010f) . Applying these minimum correction fac- 
tors to the observed column density ratios (Equation [2} leads to 
the firm lower limits listed for [C/H] m j„,[Si/H} m i n and [N/H] m j n 
abundances in columns 4-6 of Table [2] The minimum ionisation 
corrections provide firm lower limits on the abundances that do not 
depend on the ionisation uncertainties or the possibility of a multi- 
phase gas. In particular, any gas components not at an ionisation 
corresponding to IC m /„ would have the effect of raising the actual 
value of IC and thus also the actual abundance. 



Table 2. Metal Abundance and Total H Column Density 



1 The calculations in lHamann et alj 120101) apply to gas that is pho- 
toionised by a typic al quasar spectrum. We perform additional CLOUDY 
IFerland et all 19981) calculations using the inter-galactic background spec- 
trum in CLOUDY, which is based on Haardt & Madau (2005, private com- 
munication). We find that the ionisation fractions of interest in the present 
work have only negligible differences between the two calculations, e.g., 
compared to uncertainties in the measured quantities or derived ionisation 
constraints. Therefore, our analysis of the ionisation and abundances in 
J1023+5142 should apply whether the absorbers are located near the quasar 
or outside the quasar's radiative sphere of influence. 



# Zabs 


[CfflW, 


[C/H] m ,„ 


[N/H]„„„ 


[Si/H] m ,„ 


log(N(H)) cm" 2 


1 3.42865 


00.47 


>2.25 


- 


- 


17.62 


2 3.42133 


>+0.32 


>0.76 


>0.47 


- 


^17.95 


3 3.41864 


>+0.38 


>0.82 


- 


>0.34 


^17.16 


4 3.41775 


>+0.79 


>1.32 


- 


>-0.66 


^17.24 


5 3.40391 


+0.72 


^+0.02 


>0.10 


- 


18.21 


6 3.40196 


+0.50 


>0.58 


>0.32 


- 


19.13 


7 3.39889 


^+0.40 


>0.30 


>+0.01 


- 


08.19 


8 3.37985 


>+0.94 


>0.07 


>0.17 


- 


s: 18.90 


9 3.35906 


>+0.14 


^-0.86 


>+0.11 


- 


08.60 



We derive total H column densities for each NAL system from 
the H I column densities listed in column 7 of TableQ]and the best 
ionisation correction described above. We use 



log N(H) = logN(H I) - log f(H I) 



(3) 



where logN(H) is the total H column density, logN(H i) is the 
column density of H I and log f(H i) is the H I fraction used to 
obtain IC. logN(H) for each system is listed in column 7 of Table|2] 
The uncertainties in these results are dominated by uncertain- 
ties in the IC. In addition to the limited constraints provided by the 
data, a few well-studied cases have shown that individual absorbers 
can span a range of i onisations and have a range of IC values (e.g. 
lHamann et al.Ul997h ). We assume a single ionisation state for each 
absorption line system. We discuss the individual systems briefly in 



4 NOTES ON INDIVIDUAL SYSTEMS 

System 1, v = —1441 km s 

The C IV in system 1 appears line-locked with the C IV in sys- 
tem 2, as discussed further in ij l5.ll O VI is not present, or is very 
weak, implying that the ionisation is low. Further evidence for low 
ionisation is the weak C IV combined with strong H I measured in 
Ly a and Lyman y, as seen in Figure [3] Our best ionisation con- 
straint comes from an upper limit on C III, which means the best 
C/H abundance is an upper limit as well. This gas is a likely can- 
didate for host galaxy halo gas based on the weakness of the metal 
lines and the low abundances. 

System 2, v = — 1938 km s 

The C IV in system 2 appears to be line-locked with system 1, 
as mentioned above and discussed further in § 15.11 The H I could be 
shifted to a lower velocity by as much as 30 km s from the metal 
lines in this system, indicating a multi-phase gas, but heavy blend- 
ing obscures the precise shift of the lines as can be seen in Figure|4] 
The Ly a absorption line is poorly constrained. The resulting H I 
optical depth and doppler b parameter are upper limits, resulting in 
lower limits for the best estimate of C/H abundance. We constrain 
the ionisation by the relative strengths of C IV and N V, assuming 
solar abundance ratios. 

System 3, v = —2120 km s 

The H I lines in system 3 are blended with those from system 
4, but appear consistent with the metal lines, shown in Figure [5] 
Because of the relatively poor constraints on the H I absorption 
lines, the H I optical depth and doppler b parameter are upper limits, 
resulting in lower limits for the best estimate of C/H abundance. 
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System 1 



T3 
1) 




-1600 -1500 -1400 -1300 

Velocity Shift (km s~ ) 

Figure 3. Line profiles in the normalised spectrum J 1023+5 142 for system 
1. The central velocity for the Gaussian profile fit is v = —1441 km s . 
The velocity scale is with respect to the rest frame of the quasar based on 
Zem = 3.45, where negative velocities denote motion towards the observer 
and away from the quasar. The velocity range is 400 km s _1 for this and 
figures \4\ through [7] and [9] The solid curve in each panel is the Gaussian 
optical depth fit to individual lines. The dashed vertical line is the central 
velocity of the system. All of the lines used to derive or constrain column 
densities with Gaussian fits are shown in the figure. 



We constrain the ionisation by the relative strengths of C IV and 
Si IV, assuming solar abundance ratios. 

System 4, v = -2182, -2200 km s' 1 

The C IV and Si IV doublets in system 4 are fit with two 
blended Gaussian components to accommodate the asymmetric 
profile. We use the central velocity of each component to identify 
the system. The H I absorption lines are poorly constrained due to 
blending with system 3. The resulting H I optical depth and doppler 
b parameter are upper limits, resulting in lower limits for the best 
estimate of C/H abundance. We constrain the ionisation by the rel- 
ative strengths of C IV and Si IV, assuming solar abundance ratios. 
This system is broad and asymmetric, which is indicative of a wind 
or outflow feature (see § 15. It . 

System 5,v= —3121 km s 

The H I in system 5 is well constrained by Ly a . The O VI 
is strongly blended with that of system 6 as shown in Figure [6] 
The covering fraction in the doublet is ~ 0.7. The Cf = 0.7 Gaus- 
sian fits are shown as solid curves in Figure [6] System 5 appears 
to have 2 components; a narrow, optically thick component sitting 
directly on top of a broader, optically thin one. This is most clearly 



1.2 
0.8 



0.8 

0.4 
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0.4 



S °- 8 
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0.4 



0.8 
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System 2 
'C IV 1548 



C IV 1551 




4U» 



N V 1239 




A^^V'A^ 1 



N V 1243 



^^ 




-2100 -2000 -1900 -1800 

Velocity Shift (km s~ ) 

Figure 4. Line profiles in the normalised spectrum J 1023+5 142 for system 
2. The central velocity for the Gaussian profile fit is v = —1938 km s . 
The symbols and ranges are the same as in Figure|3] 



seen in Figure[6]in the longer wavelength members of the C IV and 
N V doublets, which have a much sharper central feature than their 
shorter wavelength counterparts. We constrain the ionisation by the 
relative strengths of C IV and N V, assuming solar abundance ratios. 
The partial coverage in this system indicates that it is intrinsic to the 
quasar. The partial coverage in this system and in system 6 are ex- 
amined qualitatively with the T-ratio predicted doublets, shown in 
Figure [10] and further with the point-by-point analysis, illustrated 
in Figure [TJJ Both analyses confirm similar Cf < 1 results in both 
systems (See 5 l3.2l for details). 

System 6, v = —3254 km s 

H I is well-constrained in system 6, with three mostly blend- 
free Lyman lines. The covering fraction in the doublets is Cf ~ 0.7, 
similar to system 5. The solid curves in Figure[6]represent the Cf < 
1 Gaussian fits, as for system 5. The O VI lines are blended with 
the O VI lines in system 5. The covering fraction in H I appears to 
be Cf = 1 because Ly a reaches zero intensity. We constrain the 
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Systems 4 and 3 



Systems 6 and 5 
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-2300 -2200 -2100 -2000 

Velocity Shift (km s~ ) 

Figure 5. Line profiles in the normalised spectrum J1023+5142 for sys- 
tems 3 and 4. The central velocities of the Gaussian profile fits are v = 
—2120 km s for system 3 and v = —2200 km s for system 4. The sym- 
bols and ranges are the same as in Figure[5] 



ionisation by the relative strengths of C IV and N V, assuming solar 
abundance ratios. The broad smooth shape, along with the partial 
coverage indicate that this system is part of an outflow. 

This system appears somewhat asymmetric and the f-ratio 
analysis in Figure [10] suggests further that there may be two com- 
ponents, one with partial covering near the line-centre, and a sec- 
ond broader component with complete covering in the blue wing. 
Although one component does not provide the best possible fit 
to all the lines in system 6, it is not clear that adding a second 
distinct component would provide a better characterisation of the 
actual conditions in the absorber. We test this by fitting the sys- 
tem with one and two Gaussian components, where the two com- 
ponent fit still assumes the same covering fraction in both com- 
ponents. Both fits produce similar column densities in all ions, 
A/V(C IV) = 0.15 dex, AN(N v) = 0.1 dex and AN(H i) = 0.1 dex 
in the same direction, therefore we prefer the single Gaussian fit 
in keeping with our prescription to minimise free parameters in the 




-3300 .-3200. -3100 -3000 

Velocity Shift (km s ) 

Figure 6. Line profiles in the normalised spectrum J 1023+5 142 for sys- 
tems 5 and 6. The central velocities of the Gaussian profile fits are v = 
—3121 km s for system 5 and v = —3254 km s for system 6. The solid 
curve is the Ct < 1 Gaussian fit for the C IV, N V and H I absorption lines. 
System 5 is the narrower system. The ranges are the same as in FigurefS] 



fits. Also, by using the single Gaussian fit, we ignore parts of the 
Ly a absorption which do not correspond directly to C IV absorb- 
ing gas, and therefore retain the ability to directly compare H I and 
C IV column densities for the abundance analysis. 

System 7, v = -3430,3496 km s' 1 

The H I column density is constrained as an upper limit in sys- 
tem 7 because of blending in the Lyman lines, shown in Figure [7] 



10 L. E. Simon and F. Hamann 



T3 
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System 8 



-3600 -3500 -3400 -3300 

Velocity Shift (km s~ ) 

Figure 7. Line profiles in the normalised spectrum J1023+5142 for sys- 
tem 7. This system has two blended components with central velocities of 
v = —3430 and -3496 km s . The symbols and ranges are the same as in 
Figure [3] 



resulting in lower limits for the best estimate of C/H abundance. 
We fit this broad system with two Gaussian components to better 
match the absorber shapes, and identify the system by the central 
velocities of the two components. The ionisation is constrained by 
the relative strengths of C IV and N V, assuming solar abundance 
ratios. 




-5200 -5000 -4800 -4600 -4400 -4200 

Velocity Shift (km s~ ) 

Figure 8. Line profiles in the normalised spectrum J1023+5142 for system 
8. The Gaussian profile fit has a central velocity of v = —4763 km s _I . The 
velocity range is 1100 km s _1 . The symbols and ranges are the same as in 
Figure [3] 



System 8, v = —4763 km s~ l 

The longer wavelength member of the C IV doublet in system 
8 falls on a gap between orders of the spectrograph between 6785 
and 6795 A, but the N V doublet is present in the spectrum, as is the 
shorter wavelength member of the C IV doublet. The N V doublet 
is used to determine the C/ and the Doppler b parameter for both 
doublets. This system is almost broad enough to be a mini-BAL, 
and is likely an outflow system based on the shape and strength of 
the line profile, shown in Figure[8] The H I appears to be relatively 
weak in this system compared to the metal lines, although there is 
severe blending in the Ly a forest. This blending means the H I 
absorption is poorly constrained with an upper limit, and therefore 
the best estimate for C/H abundance is a lower limit. We constrain 
the ionisation by the relative strengths of C IV and N V, assuming 
solar abundance ratios. 

We use the Gaussian fit to compare system 8 to other systems, 
but the profile of system 8 is distinctly non-Gaussian. Therefore we 
also fit the central trough of the line with a point-by-point analysis, 
shown in Figure [12] The C/H abundance found by the Gaussian 
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Figure 9. Line profiles in the normalised spectrum J1023+5142 for sys- 
tem 9. This system is a blend of three components with central velocities, 
v = —6083, -6186 and -6298 km s . Although all three components are 
fit with Gaussian profiles, only the central component is considered in the 
abundance analysis. The symbols and ranges are the same as in Figure[3] 



fit is consistent within 10% of the C/H abundance found using the 
point-by-point method. 

System 9, v = -6083, -6186, -6298 km s' 1 

System 9 has three components, but we chose to analyse only 
the central component for abundances, as the two outer components 
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Figure 10. T-predicted line profiles for systems 5 and 6. System 5 is the nar- 
rower system. The dot-dashed curve shows the smoothed C IV and N V pre- 
dicted long wavelength doublet member, based on doublet optical depth ra- 
tio from short wavelength member, assuming Cy = 1 . The actual smoothed 
data for the shorter and longer wavelength doublets are shown as the bold 
and thin solid curves, respectively. The longer wavelength data are stronger 
than the predictions, indicating partial covering, especially for the N V dou- 
blet in system 5 and in the centre of the N V doublet in system 6. 



are very poorly constrained, as shown in Figure|9] This system has 
the highest velocity shift out of the group of narrow absorption 
lines, and lies just nominally outside of the velocity shift region 
for associated lines (v ^ —5000 km s~'), at ~ —6200 km s . The 
Lyman lines could be shifted up to 20 km s from the metal lines, 
indicating a possible multiphase gas, but the line are too weak to de- 
termine their precise centroids. The weakness of the Lyman lines, 
along with blending in the Ly a forest mean the H I column densi- 
ties are upper limits, so the best estimate of the C/H abundance is 
a lower limit. We constrain the ionisation by the relative strengths 
of C IV and N V, assuming solar abundance ratios. 



5 DISCUSSION 

J1023+5142 has nine NAL systems with a range of column den- 
sities from N(H) < 10 172 to 10 191 cm" 2 , velocities from -1400 
to -6200 km s , C IV doppler b values from 7 to 150 km s , 
C IV REW(1548A) from 0.02 to 0.81 A and two systems with par- 
tial covering of either the continuum source or the broad emission 
line region (BLR), Cf ~ 0.7, which imply absorber diameters of 
^ 0.03 pc or ^ 0.8 pc (discussed below in § I5.2| |. These systems 
are generally much wea ker than those studi ed in larger statistical 
surveys of NALs, such as lVesterg aard (2003), which use lower res- 
olution data, and measure C IV REW integrated across the doublet, 
with completeness limits of 0.3-0.5 A. The NAL systems all ap- 
pear to be highly ionised; none of the systems exhibit low ionisa- 
tion species such as Si II, C II or Si III, whereas all contain C IV 
and some contain higher ionisation species such as O VI and N V. 
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Figure 11. Point-by-point covering fractions for C IV and N V in system 
6 and the centre of system 5 with step size of three resolution elements. 
System 5 is the narrower system. The solid curve is the smoothed shorter 
wavelength line, the dashed curve is the smoothed longer wavelength line, 
with their respective error spectra below. The circles represent 1 — CV at 
each step so that a point at zero flux has complete coverage, and a point at 
the continuum flux of one has no coverage. The circles are located at the 
centre of the average velocity steps. 
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Figure 12. Point-by-point covering fractions for N V in system 8 with step 
size of four resolution elements. The symbols are the same as in Figure fTTI 



Systems 5, 6, 8 and 9 exhibit high ionisation (O VI) absorption, and 
others may also have absorption at these wavelengths that is not 
observable due to blending in the Ly a forest. Systems 2-9 exhibit 
supersolar metallicities ranging from Z ^ 1 to ^ 8 Z Q . System 1 
has a slightly lower metallicity of Z ^ 0.3 Zq. We examine several 
diagnostics to estimate directly the location of each system. 

5.1 Location of the Gas 

The tight grouping and similar high metallicities (see § l5.3l for fur- 
ther discussion) for all but one (systems 2-9) of the nine C IV 



absorption line systems in J1023+5142 suggest a possible physi- 
cal connection between the absorbers. The proximity of this NAL 
complex to the quasar redshift suggests further that the physical re- 
lationship includes the quasar itself. The velocity span across the 
group is too large to be explained by a single galaxy or even a large 
cluster of galaxies. It might be consistent with some larger cosmic 
structure connected to the quasar, but then we would expect the ve- 
locity distribution to include the red side of the quasar systemic. 
A more likely explanation is that the NAL complex formed in a 
multi-component outflow from the quasar. 

There are several indirect arguments for an intrinsic origin for 
the gas in this NAL complex, i) 8 of the 9 systems have super- 
solar metallicities, discussed in detail in § 15. 3 1 below, ii) Some au- 
thors have argued th at strong O VI abs orption may indicate intrinsic 
gas near the quasar. iFox et alj ( 120081) carry out a detailed study of 
O VI absorption in 2 < z < 3 quasars and argue that logN(0 Vl) ^ 
15.0 log(cm~ 2 ) indicates an intrinsic origin, supported by evi- 
dence for partial covering in most of these systems. We measure 
N(0 Vl) in four systems in J1023+5142. Two of them (6 and 8) are 
above the intrinsic thereshold defined by Fox et al., while the other 
two (5 and 9) are very near this threshold at \ogN(0 Vl) ^ 14.5. 
iii) All systems with O VI, that is systems 5, 6, 8 and 9, also 
have strong N V compared to C IV. Strong N V, especially com- 
pared to C IV, is often (t hough not always) present in intrin- 
sic gas (e.g. Weymann et al. ( 1981); Hartquist & Snijders ( 1982); 
Hama nn et al.1 (J_997); Kuraszkiewicz & Green (2002); Fox et al. 
(2008)). iv) The presence of strong O VI and N V, especially with 
the absence of low ionisation species such as C II in these absorp- 
tion systems is consistent with gas exposed to the intense ionising 
radiation field near a quasar. 

If the NALs in J1023+5142 are intrinsic to the quasar en- 
vironment, the most likely origin is in a quasar-driven out- 
flow. Other possible intrinsic origins all have lower velocities: i) 
starburst-driven outflows typically have 100 < v < 1000 km s 
(Heckman et al. 2000), and in Seyfert galaxies have maximum 
outflow spee ds of 600 km s , and more typical speeds of 100- 
200 km s _1 JRupke et alJuOOSh . ii) other galactic/halo gas should 
have velocities near the typical velocity dispersion for such galax- 
ies (cr ~ 300 km s _1 ), iii) gas in the narrow line region of the 
quasar has typical velocities of v ^ 10 00 km s , and maxi- 
mum velocities of v < 2000 km s _1 JRuiz et al] EoOlL I2005J: 
IVeilleux et al.1 12005|). and iv) intr a-cluster galaxy motions are 
shown by Popesso & Biviano (2006) to generally have velocity dis- 



persions o" v < 1000 km s or less fo r clusters with higher num- 
bers of active galactic n uclei (Richards et al. 1999; He ckman et al.1 
l200d : IVestergaarj2003l ; lNestor et alfeOolT 

Statistical studies, iNestor et ail J2008h (see also Iwild et al.1 
J2008h ). have shown that > 43% of NALs at -750 ^ v > 
-12000 km s- 1 with REW(1548A) > 0.3 A form in high-velocity 
quasar outflows. This percentage increases to ~57% for the nar- 
rower range of —1250 ^ v ^ —6750 km s _1 , spanned by the 
NALs in J1023+5142. The percentage reaches ~72% for the nar- 
row range of —1250 ^ v ^ —3000 km s , which encompasses 
systems 1 through 4 in J102 3+5142. The s e perc entages are proba- 
bly lower fo r weaker lines ( INestor et al.1 ( 120080 and private com- 
munication). iMisawa et al.1 J2007I) also find that for C IV NALs 
with REW(1548A) > 0.056 A at velocities v < 5000 km s" 1 
the intrinsic (outflow) fraction is ^ 33% and at higher velocities, 
5000 < v < 70000 km s~\ the intrinsic fraction is > 10 - 17%. 
Nonetheless, these outflow fractions support the idea that most or 
all of the systems in this group of nine NALs in J1023+5142 form 
in a quasar outflow. 
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We search for direct signatures of quasar outflow origin via 
1) line variabi l ity, 2) partial covering and 3) br oad profiles (see 
iHamann et al.l OOlOh : lHamann & Simonl J2O10h and references 
therein for more discussion). 

1) We have only very poor constraints on the variability. We 
compare C IV and N V REW results measured from the guassian 
fits to the SDSS and Keck spectra (At rest ~ 11 months) in search 
of variability in the absorption lines. System 8 is the only individ- 
ual C IV and N V system resolved in the SDSS spectrum, while 
the weaker lines are not detected in the SDSS spectrum. System 
8 is the strongest of the nine systems, and did not vary in REW 
by more than 15% in C IV and N V between the SDSS and Keck 
observations. For the eight weaker systems, we conclude only that 
variability greater than a factor of 2 to 3 did not occur. 

2) There is partial covering in two (systems 5 and 6, Figures [6] 
andlllt out of the nine systems. Absorption lines with partial cover- 
ing of the luminosity source are attributed to gas near the quasar be- 
cause partia l covering is not expe cted to occur in intervening clouds 
or galaxies (Ha mann et alJuOlOl) . The presence of partial covering 
in these lines strongly suggests that the gas is intrinsic and located 
in the near quasar environment. 

3) The profiles of systems 4, 6, 7, 8, and 9 shown in Fig- 
ures [5] |6] |7] [8] and [9] have C IV and N V b values between 33 
and 155 km s~'and O VI b values between 42 and 150 km s . 
These b values are broad and smooth compared to the ther- 
mal widths fo r gas at the highest expec t ed temperature near 
T = 10 5 K JArnaud&Rothenflu3ll985l:lHamannetal]|l995h 
for a photoionised gas near a quasar (33 km s for H and less 
than 10 km s for C and N). They are also broader than typi- 
cal non-damped Ly a intervening C IV, N V and O VI absorp- 
tion lines, which have on average b < 12 — 14 km s for O VI 
and b < 10—12 km s~ ' for C IV and N V l Tzanavaris & Carswell 



2003 |: lBergeron & Herbert-FortfcOOSlJSchave et alj2007l ; lFox et~ 



2008). These profiles, therefore, exhibit morphologies consistent 
with formation in an outflow. 

These three characteristics, variability, partial covering, and 
broad profiles are often found together in a single object, further 
supporting the idea that each individual characteristic likely indi- 
cates an outflow. On e well studied NAL outflow in J2123-0050 
dHamann et alj 120101) is a prime example of all three, exhibit- 
ing variability, partial covering, and broad profiles that still have 
FWHMs that are as narrow or narrower than many of the systems 
inJ1023+5142. 

There is more tentative evidence for a quasar outflow origin 
in the apparent line-lock between the C IV doublets in systems 1 
and 2. Line-locking, where the difference in outflow velocities of 
two systems is exactly the velocity separation of the doublet, means 
that the lines are being radiatively accelerated directly towards the 
observer. The reality of the line-locking in this case is unclear, due 
to the difference in derived metallicities between the two systems. 
Nevertheless, the incredibly small velocity offset (§ 13. It. along with 
the v ery small probability for chance alignments (Ganguly et al. 
120031) suggest that the phenomenon may be real. The possible line- 
lock in C IV in systems 1 and 2 suggests that they are both part of 
an outflow and that these weak C IV lines play a significant role in 
radiatively driving the flow. If this is really the case here, the gas 
probably originated near the source of radiative acceleration, i.e. 
the quasar. 

Finally, we note a trend in line width with velocity shift away 
from the quasar. The narrowest lines, with FWHM ~ 20 km s 
are closest to the quasar redshift. The lines appear progressively 
broader as the velocity shift increases, with the broadest system de- 



scribed as a (narrow) mini-BAL with FWHM = 270 km s~ ' , shown 
in Figure [T] A similar phenomenon has been observed before in 
other quasars with multiple C IV absorption lines clearly forming in 
outflows dHamann etalJl 199 7; Steidel 1990; Hamann et al. 2010). 
Although it provides no direct information on the absorber loca- 
tions, the appearance of this pattern in J1023+5142 supports the 
idea that at least some of the systems form in a quasar outflow. The 
tight grouping of all nine of the sys tems also suggests a relationship 
between them. lGangulv et alj d2003l) determine that the probability 
of six similarly grouped NALs in the quasar RX Jl 230.8+01 15 all 
forming in intervening (uncorrelated) gas at similar velocity shifts 
is extremely small. The similarities between those NALs and the 
NAL complex in J1023+1542 implies a similarly small probability 
for all nine NALs in this complex forming independently in inter- 
vening gas. Although there could be up to several interlopers in the 
NAL complex of J1023+5142 that might form in nearby galaxies 
in the line of sight, the density of these galaxies required to form 
all of the absorbers in the complex is beyond any expectations of 
cluster density at this redshift. 

Overall, we conclude that at least six out of the nine systems 
originate in a highly structured outflow driven by the quasar, be- 
cause they exhibit one or more of the following properties: par- 
tial covering, broad profile shapes, large line strengths, tight group- 
ing with other systems, and proximity to the quasar redshift. Sys- 
tems 5 and 6 are the most likely outflow candidates because they 
exhibit partial covering as well as several of the other properties 
listed above. Systems 4 and 8 are likely outflows because of their 
strong, broad, asymmetric and smooth profiles and systems 7 and 9 
are probably outflows because of their broad and smooth shapes. 
Systems 1 through 3 are more ambiguous in origin, with nar- 
row widths, complete covering, lower velocity shifts and smaller 
strengths. However systems 1 and 2 exhibit line-locking, which 
could be evidence of an outflow. Ultimately, we find strong evi- 
dence that systems 4, 5, 6, 7, 8 and 9 are part of a quasar outflow, 
whereas systems 1, 2, and 3 could consist of intervening gas from 
the IGM or other galaxies in the line of sight. 



5.2 Outflow Properties 

As described in § 15.11 the evidence suggests that the majority of 
absorption lines in this grouping are part of a complex quasar out- 
flow. This flow must be highly structured, with at least 6 and as 
many as 9 distinct absorbing structures along the line of sight. The 
velocities in the 6 most secure outflow systems range from -2120 
to -4760 km s _1 . Several of these systems (4, 6, 7 and 8) also have 
super-thermal line widths, indicative of large turbulence or strong 
radial velocity sheer across the outflow structure. 

Some of the outflow structures, represented by systems 5 and 
6, must be spatially small to produce partial covering of the back- 
ground emission source. These lines lie on top of the very weak 
C IV broad emission line. Therefore, nearly all (>80%) of the flux 
beneath these lines is continuum emission and any partial covering 
below Cf = 0.80 can be ascribed to the continuum source and not 
the much larger BLR. We estimate the diameter of the accretion 
disk continuum source at 1550 A to be d ~ 0.03 pc and the diam- 
eter of the C IV broad line region to be d ~ 0.8 pc, based on the 
scaling relational in lHamann & Simonl dJ2010f> . To partially cover 



2 We estimate the luminosity from the rest-frame flux at 1450 A measured 
in the SDSS spectrum, which, combined with the luminosity distance and 
the bolometric correction factor L=3.4*vL v (1450), gives AL A (1450A). We 
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the emission source, the absorbing clouds should have characteris- 
tic sizes similar to or less than the BLR diameter, and possibly even 
less than the accretion disk diameter. 

If the absorbers are discrete clouds, their small sizes and sub- 
stantial velocity dispersions should lead to fa irly rapid dissipa- 
tion in the absence of an external pressure (see lHamann & Simonl 
02010]) ). In particular, the characteristic size of d ~ 0.03 pc or 
d ~ 0.8 pc combined with b = 45 km s in the partial covering sys- 
tem 6 indicate a dissipation time of roughly f,/, s ~ d/b ~ 660 yrs or 
tdis ~ 17,400 yrs. At the measured velocity of v = —3254 km s , 
this gas component would travel just ~ 2 pc or ~ 60 pc before dissi- 
pating. A thorough discussion of the creation and survival of these 
absorbing structures is beyond the scope of this paper. However, 
these simple arguments suggest that at least some of the outflow 
components we measure are very near their point of creation. 

It is useful to compare the basic properties of this NAL out- 
flow to BALs. The ionisations in both types of outflows are simi- 
lar, with very little low ionisation gas (e.g. C III). The outflow ve- 
locities of the NALs in J1023+5142 are lower, < 6200 km s -1 , 
than the typical outflow velocities in BAL s, which can reach 
up to > 20,000 km s _1 JKorista et allll993h . but they do over- 
lap. The NALs have total H column densities (N(H) < 10 172 to 
10 cm , individual values listed in Tabled, more than 1000 
times lower than typical BAL H column densities (N(H) ^ 10 20 — 
10 cm , and probably higher). By definition, these NALs also 
have FWHMs around 1000 times narrower than typical BALs, and 
much smaller REWs as well. Nonetheless, NALs like this might be 
part of the same general outflow phenomenon as BALs, viewed at 
different angles dElvisfeOOOklGangulv et ai1l200lh . 

This complex of weak NAL outflows appears to be dramati- 
cally different from typical BAL outflows, and constitutes a nearly 
unexplored part of the quasar outflow phenomenon, with a range 
of physical parameters and kinematics more complex and varied 
than previously thought. It is well-known that NALs are a common 
feature of quasar spectra. Previous surveys have found that 40% of 
quasars have C IV NALs, in particular 25% h ave strong C IV NALs 
within v > —5000 km s _1 dVestergaardll2003l) , 60% have quasar- 
driven outfl ows in some form, either BALs, NALs, or something 
in-be tween (Ganguly & Brotherton 2008; R odriguez Hidalgo et al.l 
UOIOI) . and includ ing high-velocity outflows raises the percentage 
to 70% flMisawa et alj|2007h . If the coverage fraction of these out- 
flows is less than 100%, which is likely, they could be ubiquitous 
in the near quasar environment, and could potentially play an influ- 
ential role in the physical processes occurring therein. 

Finally, we would like to understand what role the NAL out- 
flow in J1023+5142 might have in feedback to galaxy evolution. 
The low speeds and small column densities, e.g., compared to BAL 
flows, suggest that its feedback contribution is negligible. How- 
ever, there are large uncertainties relating to the outflow location 
and geometry. For one particular NAL o utflow at a derived radial 
distance of ~ 5 pc, lHamann et al.l ( 120101) estimate that the kinetic 
energy yield is several orders of magnitude smaller tha n that nec- 
essary to influence feedback. At the opposite extreme, iMoe et al.l 
(2009) argue that the feedback contribution is significant for an- 
other NAL outflow at a derived distance of ~ 2 — 5 kpc. The loca- 



measure the C IV emission line FHW M in the SDSS spectrum, and using 
Equation 7 in lVestergaard & Peterson! (200g), derive a black hole mass of 
\ogMgH = 9.8 Mi. Based on these values, we calculate an Eddington lu- 
minosity fraction of 0.8. The black hole mass and Eddington luminosity 
fraction are then used in the scaling relations in Hamann & Simon (2010) 
to calculate the size of the continuum and broad line emission regions. 



tion of the NAL outflow in J1023+5142 is not known well enough 
to make these estimates. A more sensitive search for variability in 
these NALs could be very h elpful for refin ing both the location 
and the total energy yield (see lHamann et al . ( 2010) and references 
therein). 



5.3 Metallicity 

We find greater than or consistent with solar abundances in all 
of the absorption systems in J1023+5142 except in system 1, in 
agreement with p revious studies of narrow associated absorption 
at lower redshifts (Petitjean & Srianand 1999; Ha mann et al. 2001; 



iD'Odorico et alj |2004 iGabel et al.l 120061) . These high metallic! 
ties are consistent with the re sults of other studies of intrinsic gas 
as well, including BEL gas ( lHamann et al.l 120020 and therefore 
consistent with our interpretation that the gas is intrinsic to the 
quasar. Intervening absorbers generally have very low metallici- 
ties, with Z no more than a few hundredths solar, although there 
are cases where high metallicity intervening gas has been observed 
(Prochaska et al. 2006; Schave et al. 2007). We argue that the high 
metallicities found in 8 of the 9 systems in this quasar are consistent 
with locations near the quasar, however, we do not rely solely on 
this argument to determine the gas location. Instead, we consider 
that high metallicities could be a general phenomenon found in all 
gas in the quasar host environment (Prochaska & Hennawi 2009). 
The high metallicities of the NAL systems in J1023+5142 
require that its host galaxy had vigorous star formation in the 
epoch before the quasa r was observable, le ading to metal-rich gas 
in the quasar outflows JFalomo et al.l l2008). This evidence, along 
with previous studies of BELs lead us to conclude that the gener- 
ally accepted paradigm of quasar-host galaxy evolution is correct, 
where a major merger leads to a vigorous burst of star formation, 
which then funnels gas to the centre of the galaxy and ignites a 
quasar that eventually blows out obscuring gas and dust to become 
visibly luminous (Perez-Go nzalez et al . 2008; Ho pkins et al . 2008; 
iRamos Almeida et alj 2009). However, larger samples are needed 
to examine the full range of NAL properties and study their rela- 
tionships to quasar outflows and host galaxy environments. Mea- 
surements at high redshifts are particularly valuable because this is 
the main epoch of host/massive galaxy formation when the NAL 
gas might have a close relationship to ongoing or recent star forma- 
tion in the hosts. 



6 SUMMARY 

We use NALs to improve our understanding of the evolutionary 
relationship between the central black hole and its host galaxy 
through the study of their location, origin and abundance informa- 
tion in high redshift quasars. Here, we examine the properties of 
nine NAL systems in the quasar J1023+5142 and find N(H) from 
< 10 17 - 2 to 10 19 ' 1 cm" 2 , velocities from -1400 to -6200 km s" 1 , 
C IV doppler b values from 7 to 150 km s~\ C IV REW from 0.02 
to 0.81 A and two systems with partial covering of either the con- 
tinuum or the BLR, at the level of Cf ~ 0.7, which imply absorber 
diameters of ^ 0.03 pc or ^ 0.8 pc. 

The NAL systems all appear to be highly ionised; none of the 
systems exhibit low ionisation species such as Si II, C II or Si III, 
whereas all contain C IV and some contain higher ionisation species 
such as O VI and N V. 

The C IV absorption NALs are tightly grouped, suggesting that 
they have a physical relationship to one-another, and the proximity 
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of the NAL complex to the quasar redshift suggests that the phys- 
ical relationship includes the quasar itself. The range in velocity 
across the complex is larger than can be easily explained by a sin- 
gle galaxy or even by a large cluster of galaxies. A more likely 
explanation is that the NAL complex formed in a multi-component 
quasar-driven outflow. 

We estimate directly the location of each system in 
J1023+5142 through the use of several diagnostics and find strong 
evidence (partial covering, broad and smooth profiles compared to 
thermal widths, velocities greater than galaxy dispersion velocities, 
supersolar metallicities) that systems 4, 5, 6, 7, 8 and 9 are part of 
a quasar outflow. Systems 1, 2, and 3 have more ambiguous ori- 
gins because they exhibit narrow widths, lower velocity shifts, and 
system 1 has a lower metallicity, so these systems could consist of 
intervening gas from the IGM or other galaxies in the line of sight. 

Systems 2-9 in J1023+5142 exhibit supersolar metallicities 
ranging from Z ^ 1 to ^ 8 Zq. System 1 has a lower metallicity of 
Z ^ 0.3 Zq. The high metallicities are consistent with scenarios of 
galaxy and black hole formation and evolution. 

The NALs in outflows appear to be part of a related outflow 
complex, which is very different than other known outflow regions 
such as BAL outflows, and constitutes a relatively unexplored part 
of the quasar outflow phenomenon. The outflows in J1023+5142 
could be important for feedback between the black hole and the 
host galaxy, depending on the radial distance of the gas from the 
quasar. 

The narrow widths of NALs mean that detailed studies of in- 
dividual objects like this are the only way to make progress in un- 
derstanding this type of outflow. Variability studies could be useful 
to add more examples of NAL outflows to the current sample avail- 
able for similar detailed analysis. 

We will add significantly to this sample in future work, in- 
cluding detailed studies of the full range of NAL properties in 24 
quasars at high redshift, during the main epoch of host/massive 
galaxy formation. 
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